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1.  Introduction 

SASPRO,  an  acronym  standing  for  5pare  and  Server  Provisioning, 
is  a versatile  FORTRAN  package  that  gives  provisioning  levels  of  spares 
inventory  and  repair  capacity  required  to  support  a population 
of  randomly  failing  items  which,  upon  failure,  are  (1)  dispatched  to 
the  repair  facility,  and  (2)  replaced  by  a spare  if  one  is  available. 
This  paper  describes  in  detail  the  problem  environment,  the  program  op- 
tions, the  input  required  to  run  the  program,  and  the  output  provided 
by  the  program.  Sample  runs  are  also  shown  for  each  of  the  program 
options. 

2.  Problem  Environment 

A population  of  items  containing  certain  key  parts  (for  example, 
a fleet  of  aircraft,  a fleet  of  ships,  a group  of  machines)  randomly 
falls  and  requires  repair.  Spare  parts  are  also  needed  so  that  upon 
failure,  the  spare  can  be  utilized  to  replace  the  failed  part  and  the 
item  put  back  into  service.  It  is  desired  to  determine  how  many  spares 
and  how  many  repair  channels  are  required  to  support  the  system  at  a 
desired  service  level  while  minimizing  costs. 
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The  system  is  shown  schematically  in  Figure  1.  We  consider  only 
a single  part-type  at  a time,  which  has  its  own  spares  pool  and  dedicated 
repair  channels.  For  example,  for  a fleet  of  gas  turbine  propelled  ships, 
the  gas  turbine  engine  has  two  components — a gas  generator  and  a power 
turbine.  Each  must  have  dedicated  repair  channels  and  its  own  spares 
pool.  Thus  SASPRO  would  treat  each  component  in  turn,  being  utilized  to 
provision  first  for  a population  of  gas  generators  and  then  for  a popula- 
tion of  power  turbines. 

When  a unit  in  the  operating  population  fails,  a spare  is  requested 
at  the  same  time  the  unit  is  dispatched  for  repair.  If  a spare  is  not 
available,  the  request  is  backlogged  and  units  coming  out  of  repair  are 
used  in  removing  the  backlog.  When  there  is  no  backlog  of  requests  for 
spares,  units  coming  out  of  repair  go  into  the  spares  inventory.  Repair 
times  as  well  as  failure  times  are  treated  as  random  variables  and  with 
the  proper  assumptions  (to  be  mentioned  below) , this  stochastic  process 
can  be  readily  modeled  as  a finite  source  queueing  system,  often  referred 
to  as  "the  machine  repair  problem  with  spares,"  which,  in  addition,  also 
fits  a two-stage  cyclic  queueing  model.  Thus  SASPRO  uses  a standard 
queueing  model  for  the  stochastic  process  [see  GROSS,  KAHN,  and  MARSH 
(1977)]. 

The  assumptions  required  for  using  SASPRO  are  that  times  to  fail- 
ure and  repair  times  are  exponentially  distributed  random  variables. 

These  assumptions  allow  the  employment  of  the  standard  finite  source 
queueing  theory  to  determine  probabilities  of  various  numbers  of  units 
in  repair  at  any  given  time.  From  this,  system  service  levels  (number 
of  units  operating,  availability  of  spares,  etc.)  can  readily  be  computed. 
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In  order  to  achieve  a specified  service  level,  certain  combinations  of  spares 
and  servers  (repair  channels)  are  required.  Using  costs  associated  with 
purchasing  and  holding  spares,  and  costs  associated  with  building  and  op- 
erating repair  channels,  SASPRO,  through  a heuristic  optimization  routine, 
finds  the  ’’best"  combination  of  spares  and  servers  that  meets  the  service 
level  constraint.  [See  Gross,  et  at.  (1977)  for  specific  details.] 

3.  Modes  of  Operation 

The  program  has  two  modes  of  operation,  dynamic  and  static.  The 
former  is  advised  for  initial  year  provisioning  when  population  sizes, 
failure  rates,  and  repair  times  may  be  changing  significantly.  Population 
size  changes  may  occur  because  units  are  put  into  operation  gradually, 
thus  building  up  to  a full  strength  population  over  a period  of  several 
years.  For  example,  it  was  anticipated  to  build  a fleet  of  256  gas  tur- 
bine powered  ships,  starting  in  the  first  year  with  ten  ships  and  building 
up  to  full  strength  over  a ten-year  period.  Because  of  new  technology, 
engines  on  ships  introduced  in  the  later  years  were  expected  to  have 
smaller  failure  rates  (be  more  reliable)  while  due  to  learning,  repair 
times  were  also  expected  to  be  smaller  in  the  later  years. 

In  designing  support  systems  for  which  it  is  necessary  to  deter- 
mine the  number  and  location  of  depots  the  static  mode  is  useful  [see 
GROSS  and  PINKUS  (1978) ] . In  this  situation  the  population  is  at  full 
strength,  technological  advances  and  learning  are  complete,  and  condi- 
tions are  very  close  to  static.  Running  times  and  input  requirements 
are  greatly  reduced  when  operating  in  the  static  mode. 

The  dynamic  mode  allows  for  changing  population  sizes,  failure 
rates,  and  repair  times,  as  well  as  for  changing  costs  on  a year  by  year 
basis.  A set  of  input  information  is  required  for  each  year  in  the  plan- 
ning horizon.  An  item  repaired  and  placed  back  into  spares  inventory  (or 
operation)  in  year  i is  assumed  to  have  the  same  failure  rate  as  a new 
item  Introduced  in  year  1 . 

All  costs  are  turned  into  an  equivalent  beginning  of  year  payment 
and  then  discounted  according  to  where  the  year  is  in  the  planning  horizon 
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so  that  at  any  year  i the  program  gives  the  present  worth  of  the  cumu- 
lative sum  of  the  discounted  costs  up  to  and  including  year  i . The 
final  value  for  the  last  year  is  then  the  present  worth  of  the  sum  of 
discounted  costs  over  the  entire  planning  horizon. 

The  costs  that  are  considered  in  SASPRO  are  purcnase  costs,  sal- 
vage values,  and  annual  operating  costs  associated  with  spares  and  repair 
channels;  unit  transportation  and  repair  costs;  and  component  improvement 
investment  costs.  Purchase  costs  and  salvage  values  are  in  dollars  per 
spare  or  repair  channel,  and  salvage  values  are  not  realizable  until  the 
end  of  the  planning  horizon,  even  though  spares  or  repair  channels  are 
retired  prior  to  that.  Operating  costs  of  each  channel  and  holding  costs 
for  spares  are  in  dollars  per  year  per  spare  or  repair  channel  and  any 
spare  or  channel  purchased  in  year  1 is  assumed  to  incur  annual  costs 
until  the  end  of  the  planning  horizon,  even  if  retired  earlier.  Transpor- 
tation and  repair  costs  are  in  dollars  per  unit  per  repaired  item  and  com- 
ponent improvement  program  (CIP)  costs  are  in  dollars  per  year. 

The  assumptions  that  salvage  values  are  not  received  until  the  end 
of  the  planning  horizon  and  that  operating  or  holding  costs  are  not  re- 
duced when  spares  or  repair  channels  are  retired  early  are  necessitated 
by  the  heuristic  optimization  algorithm  employed.  Further,  if  a spare  or 
repair  channel  is  retired  during  the  planning  horizon  and  is  required 
again  a few  years  later,  it  must  be  repurchased.  Again,  this  assumption 
is  required  because  of  the  heuristic  optimization  algorithm. 

In  the  static  mode  no  such  problems  arise  since  spares  and  channels 
are  not  added  or  retired.  Costs  for  this  mode  of  operation  are  converted 
to  expected  equivalent  end-of-year  payments  over  the  planning  horizon. 

4.  Service  Level  Constraint  Options 

There  are  two  options  available  for  specifying  service  performance. 

The  first,  referred  to  as  spares  availability,  sets  a limit  on  the  percentage 
of  requests  for  spares  that  are  met  from  on-shelf  spares  inventory  (also 
called  fill  rate);  that  is, 
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Nunber  of  Requests  for  Spares  per  Year  Honored  Inmedlately 
Number  of  Requests  for  Spares  per  Year 


where 


is  specified  by  the  user  and  0 < A < 1 . 


The  second  criterion  for  service  performance,  called  "fleet"  avail- 
ability, sets  a level  for  the  percentage  of  time  a certain  portion  of  the 
population  desired  to  be  in  operation  is  actually  operating;  that  is, 

Pr{^  SM  units  are  up}  ^ A . (2) 

Both  6 and  A are  specified  by  the  user  when  0 < 3 ^ 1 , and  M is 
the  population  (fleet)  size  excluding  spares. 


Suppose,  for  example,  we  wish  to  have  100  machines  in  operation 
(M“100).  When  a machine  fails,  a spare  machine  is  "plugged  in"  if  one 
is  available.  We  might  specify  that  a service  level  constraint  be  (1) 
the  percentage  of  requests  for  spare  machines  filled  immediately  from 
on-hand  spares  is  at  least  90%  (option  1:  A =.9),  or  (2)  at  least  95% 
of  the  machines  are  operating  85%  of  the  time  (option  2:  B=.95,  A=.85) 


5.  Input  Data 

Table  I shows  the  data  that  are  required  as  input  for  SASPRO. 


Most  input  parameters  are  self-explanatory  but  a few  require  further 
comment.  The  A shown  in  Equations  (1)  and  (2)  is  AVL,  while  BETA  is  the 
3 shown  in  Equation  (2). 


The  initial  values  CO  and  YO  must  be  read  in  for  number  of  ser- 
vers and  spares,  respectively.  In  the  dynamic  mode,  after  the  first  year 
I the  program  uses  the  previous  year' s values  for  CO  and  YO  as  initial 

values,  even  though  a CO  and  YO  must  be  specified  on  the  card  set  for 
every  year  in  the  planning  horizon.  The  closer  the  initial  values  are  to 
the  final  values  (determined  by  SASPRO),  the  fewer  iterations  of  the  heu- 
ristic optimization  algorithm  are  required.  However,  one  may  use  CO  • 

* YO  ■ 1 if  so  desired. 

I 
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TABLE  I 

INPUT  REQUIREMENTS 


Variable 

Name 

Description 

Symbol 

on 

Printout 

AVL 

Desired  Availability 

AVL 

BETA 

Desired  Percent  of  Population  Up 

BETA 

C 

Initial  Value — Number  of  Repair  Channels 

CO 

CIC 

Carrying  Cost  per  Spare  ($/yr/spare) 

CIC 

CIPC 

Component  Improvement  Cost  ($/yr) 

CIPC 

CPSER 

Repair  Channel  (Server)  Purchase  Cost  ($/channel) 

CPSER 

CPSP 

Spare  Purchase  Cost  ($/spare) 

CPSP 

H 

Operating  Hours  per  Year  per  Item  (hrs) 

H 

KEYWD 

Mode  Option  Indicator 

^ Static 

KEYWD 

KZ 

Service  Criterion  Option  Indicator  1^?'* 

(=1:  Fleet 

Avl 

Avl 

KZ 

NYEARS 

Planning  Horizon  Length  (yrs) 

YEARS 

OCSPSER 

Operating  Cost  of  a Channel  ($/yr /channel) 

OCPSER 

R 

Yearly  Interest  Rate 

RATE 

RM 

Population  Size 

M 

RMTBR 

Mean  Time  Between  Removals  (hrs) 

MTBR 

ST 

Average  Turn  Around  Time  (days) 

1/MU 

SVSER 

Salvage  Value  of  a Channel  ($/channel) 

SVSER 

SVSP 

Salvage  Value  of  a Spare  ($/spare) 

SVSP 

URC 

Unit  Repair  Cost  ($/unit) 

URC 

UTC 

Unit  Transportation  Cost  ($/unlt) 

UTC 

Y 

Initial  Value — Number  of  Spares 

YO 

A set  of  cost  Inputs  (CIC,  CIPC,  CPSER,  CPSP,  OCSPSER,  SVSER, 

SVSP,  URC,  UTC)  is  required  for  each  year  in  the  horizon  in  the  dynamic 
mode.  This  allows  one  to  account  for  Inflation  and  technological  inno- 
vations. The  Component  Improvement  Program  Cost  (CIPC)  is  the  annual  ex- 
penditure required  to  achieve  a given  MTBR  schedule  (the  MTBR  which  must 
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be  Inputted  for  each  year  in  the  horizon)  for  the  dynamic  mode  of  opera- 
tion, or  to  maintain  the  MTBR  achieved  when  operating  in  the  static  mode. 

The  MTBR  value  is  the  actual  mean  time  to  failure  of  each  unit 
when  operating  continuously  If  items  do  not  operate  continuously  but 
are  required  for,  say,  only  H hours  per  year  on  the  average,  the  mean 
failure  rate  actually  used  in  the  queueing  model  portion  of  SASPRO  is  low- 
ered accordingly.  If  items  do  operate  around  the  clock,  H = 365  x 24  = 
8760.  If,  for  example,  a population  of  items  has  an  MTBR  of  1000  hours 
but  is  called  upon  to  operate,  on  the  average,  only  half  the  time  (H  = 

4380  hours),  the  effective  MTBR  used  in  the  program  is  raised  to  2000 
hours  (failure  rate  cut  in  half).  The  user  specifies  H and  MTBR  and 
SASPRO  automatically  makes  the  adjustment.  The  reader  is  referred  to 
BARZILY,  GROSS,  and  KAHN  (1977)  for  a discussion  of  the  adequacy  of  this 
procedure  to  account  for  noncontinuous  operation.  The  above  reference 
also  discusses  the  SASPRO  assumptions,  when  operating  in  the  dynamic  mode, 
that  (1)  the  population  attains  instantaneous  steady-state  each  year  at 
average  values,  and  (2)  the  population  consists  of  nonidentical  units 
(with  respect  to  mean  failure  and  repair  time),  which  are  treated  as  iden- 
tical by  weighted  averaging. 

The  parameters  KEYWD  and  KZ  are  the  option  flags.  Setting  KEYWD=1 
puts  SASPRO  in  the  dynamic  mode;  otherwise  it  operates  in  the  static  mode. 
Designating  KZ=1  sets  the  service  level  on  fleet  availability;  otherwise 
the  service  level  constraint  is  on  spares  availability.  A 3 must  be 
specified  when  KZ=1  . 

Table  II  and  Figure  2 give  the  card  layout  required  for  the  input 
information.  There  are  eight  cards  needed  for  the  static  mode  and  six 
plus  two  cards  for  each  year  in  the  planning  horizon  required  for  dynamic 
mode  operation.  The  input  requirements  for  static  mode  operation  are 
similar  to  those  required  for  a one-year  planning  horizon  dynamic  run. 
However,  the  output  cost  values  given  in  the  static  mode  are  the  expected 
end  of  year  payments  adjusted  over  an  NYEAR  life,  while  the  costs  of  a 
• single  year  dynamic  mode  run  are  the  present  worth  of  expenditures  for 

that  year. 
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CARD  LAYOUT  FOR  INPUT 


Card 

Number 

Input  Data  Parameter (s) 

Format 

Columns 

1 

Title  (any  desired  by  user) 

— 

1-80 

2 

NYEARS 

12 

1-2 

3 

R 

F8.5 

1-8 

4^ 

AVL,  BETA 

F8.5,  F8.5 

1-8,  9-16 

5^ 

KZ 

12 

1-2 

6" 

KEYWD 

12 

1-2 

See  Figure  2:  One  set  required  for 
each  year  in  dynamic  planning  hori- 

8 

zon;  one  set  only  for  static  mode. 

^For  Spares  Avail  Option, 


BETA  may  be  set  at  any  value. 


b 


KZ 


0 -»■  Spares  Availability 

1 -*■  Fleet  Availability 


^KEYWD 


0 ■>  Static  Mode 

1 -►  Dynamic  Mode 


6.  Output  from  SASPRO 

SASPRO  gives  the  heuristic  optimum  combination  of  spares 
and  repair  channels  needed  to  meet  the  service  level  constraint  and  also 
provides  the  costs  associated  with  this  solution.  For  the  static  opera- 
tion mode  there  is  a single  line  of  output  with  the  total  costs  shown 
being  the  expected  equivalent  annual  expenditure  over  the  NYEAR  system 
life.  For  the  dynamic  mode  of  operation  there  is  a line  of  output  for 
each  year,  the  cost  outputted  being  the  expected  present  worth  of  the 
cumulative  sum  of  discounted  costs  up  to  and  including  that  year.  Also 
given  as  output  are  the  heuristic  optimum  combinations  of  servers  and  spares 
the  average  system  failure  rate,  which  in  the  static  mode  is  the  same  as  the 
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inputted  failure  rate  calculated  from  the  MTBR  and  H values  read  in,  and 
in  the  dynamic  mode  is  a weighted  average  of  the  units  in  the  population 
which  were  introduced  and  repaired  in  various  years  at  different  values;^ 
the  average  number  of  units  repaired;  and  the  actual  availability  achieved 
(always  ^ AVL) . Another  output  quantity  shown  is  ASTAR,  the  percentage 
of  time  the  population  is  called  upon  to  operate  (ASTAR  = H/8760) ; this 
serves  as  a check  on  the  H value  put  in.  The  output  quantities  with 
definitions  are  shown  in  Table  III. 


TABLE  III 
OUTPUT  QUANTITIES 


Name 

Description 

YR 

Actual  year  represented 

M 

Population  size  year  i (from  input) 

FR 

Failure  rate  of  a unit  purchased  or  repaired  in  year  i 
(failures/day  = 

FRBAR 

Average  failure  rate  of  a typical  unit  (failures/day  = 
weighted  average  of  various  units  purchased  or  repaired 
in  all  years  up  to  and  including  i) 

ASTAR 

Average  percent  of  time  population  is  called  upon  to 
operate  (H/8760) 

C 

Heuristic  optimum  number  of  repair  channels  required  in 
year  i 

Y 

Heuristic  optimum  number  of  spares  required  in  year  i 

RBAR 

Average  number  of  units  repaired  in  year  i 

COST 

Costs  expended  in  year  i dynamic  mode  or  equivalent 
yearly  average  expenditure  in  static  mode 

PRES  WORTH 

Present  worth  of  sum  of  discounted  costs  up  to  and  in- 
cluding year  1,  dynamic  mode;  same  as  COST  for  static 
mode 

1 


See  Gross,  Kahn,  and  Marsh  (1977). 
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7.  Sample  Runs 

We  illustrate  the  model  by  showing  sample  runs  for  each  mode  (dy- 
namic and  static)  for  each  service  level  option  (spares  availability  and 
fleet  availability).  The  first  run  is  a dynamic,  fleet  availability,  20- 
year  planning  horizon  case.  Following  this,  a similar  dynamic  run  is  per- 
formed for  spares  availability.  Next  a static,  20-year  system  life  spares 
availability  situation  is  considered.  Finally,  a static,  20-year  system 
life  fleet  availability  case  is  run. 

A listing  of  the  input  cards  for  these  runs  is  given  in  Figure  3. 
For  the  four  cases,  there  are  a total  of  108  data  input  cards  (2 [6  + 

2 X 20]  + 2[8])  . 

The  associated  output  for  the  cases  is  given  in  Figure  4. 

8.  Intermediate  Output  and  the  Heuristic 
Optimization  Algorithm 

Also  provided  as  output  are  intermediate  values  of  Y and  C 
which  "step  up"  from  YO  and  CO  showing  the  operation  of  the  algorithm 
at  each  iteration.  Briefly,  the  algorithm  works  as  follows. 

First,  a and  are  computed,  where  is  a function  of 

the  purchase  cost,  operating  cost,  and  salvage  value  of  a repair  channel 
and  is  a function  of  the  purchase  cost,  carrying  cost,  and  salvage 

value  of  a spare.  In  the  dynamic  mode,  the  functions  for  and  C2  , 

year  i,  are 

C = CPSER  + OCSPER  " SVSER  [ , ,,1  , 

L R(l+R)  J L(1+R)*^“^ 

(3) 

S - --  ^ -C  - -P  [-^]  . 

where  K is  the  total  number  of  periods  in  the  dynamic  planning  horizon, 

1 the  current  year,  R the  Interest  rate,  and  the  costs  are  as  defined 
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Figure  3.  Input  card  listing  for  sample  runs. 
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in  Table  I.  The  first  bracket  term  brings  the  annual  costs,  OCSPER  and 
CIC,  to  a beginning  of  year  i equivalent  cost,  while  the  second  bracket 
term  brings  the  salvage  value  to  a beginning  of  year  i equivalent  term; 
that  is,  the  bracket  terms  are  the  present  worth  factors  for  a series  pay- 
ment and  end  of  horizon  payment,  respectively.  Note  it  is  assumed  that 
if  a spare  or  repair  channel  is  purchased  in  year  i,  the  annual  costs  are 
Incurred  through  the  end  of  the  planning  horizon  and  salvage  cannot  be 
obtained  until  the  end  of  the  horizon,  even  if  removal  occurs  sooner. 

The  algorithm  forms  a ratio  (call  A ) of  ^^^^2  ’ 

depending  on  the  relative  magnitudes  in  such  a way  that  the  ratio  is  > 1 • 
Then  given  a pair  of  values  C,Y  (to  start  year  i,  and  ^ are 

used)  the  availability  is  computed.  If  it  is  below  the  desired  level  and 
if  A = C^/C2  , then  for  an  equal  dollar  expenditure  A repair  channels 

or  one  spare  can  be  added.  Availability  is  calculated  for  both  cases 
(adding  A repair  channels  or  one  spare)  and  the  case  yielding  the  high- 
er availability  becomes  the  new  C,Y  pair.  The  algorithm  repeats  until 
the  desired  availability  is  met.  Upon  exceeding  the  desired  availability, 
a backoff  procedure  is  utilized.  If  feasibility  was  reached  by  adding  A 
channels,  the  algorithm  first  attempts  to  remove  a spare  and  then  channels 
are  removed  one  at  a time  to  see  if  a cheaper  solution  exists  near  the 
boundary.  If  feasibility  was  reached  by  adding  a spare,  again  one-at-a-tlme 
removal  of  channels  is  tried.  Had  A = words  channel  and  spare 

would  be  reversed  in  describing  the  algorithm. 

When  the  initial  values  of  C and  Y for  year  i exceed  the 
availability  desired,  the  algorithm  immediately  goes  into  a backoff  mode, 
trying  to  remove  spares  and  channels  one  at  a time,  starting  with  the  more 
expensive  (larger  value)  first. 

The  algorithm  uses  only  and  C2  • The  other  costs  are  used 

in  calculating  the  total  year  i expenditures  even  though  they  are  not 
used  in  the  algorithm.  The  costs  URC,  UTC,  and  CIPC  are  assumed  year 
end  expenditures  and  are  multiplied  by  1/ (1+R)  to  bring  them  back  to 
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a year  beginning  cost.  The  URC  and  UTC  costs  are  multiplied  by  R and 
then  added  to  C^x(nuinber  of  new  channels  added)  + C^’^Cnumber  of  new 

spares  added)  + CIPC  to  yield  the  total  beginning  of  year  cost  for  year 
i (call  ) which  is  outputted.  Letting  a = 1/(1+R)  , the  discount- 

ing factor,  the  final  cost  figure  given  in  the  output  for  year  i is  the 
present  worth  of  the  cumulative  sum  of  equivalent  beginning  of  year  costs, 


up  to  and  including  year  i;  namely. 


In  the  static  mode  the  algorithm  works  in  the  same  way,  except 
the  functions  and  are  changed  to  reflect  all  costs  as  average 

expenditures  per  year.  Thus,  the  purchase  costs  and  salvage  values  are 
multiplied  by  sinking  fund  and  capital  recovery  factors,  namely. 


In  this  mode  total  equivalent  expenditures  per  year  are  calculated 
simply  by  adding  (URC  + UTC)R  and  CIPC  to  Cj^x(c)  + . This  is 

then  the  value  which  shows  as  both  COST  and  PRES  WORTH  on  the  output,  the 
PRES  WORTH  column  being  meaningless  in  the  static  mode. 

A sample  of  intermediate  output  is  shown  in  Figure  5 for  the 
first  year  of  sample  run  1.  Shown  are  the  failure  rate  for  year  i (RLAM) , 
average  population  failure  rate  for  year  1 (AMTBR) , average  turn  around 
(repair)  time  (ST),  availability  for  the  particular  combination  of  C and  Y , 
average  queue  size  at  repair  depot  (LQ) , and  average  number  of  units  in 
repair  (L) . 
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To  cope  with  the  expanding  technology,  our  society  must 
be  assured  of  a continuing  supply  of  rigorously  trained 
» and  educated  engineers.  The  School  of  Engineering  and 

Applied  Science  is  completely  committed  to  this  ob- 
jective, 
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